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ABSTRACT 

There is evidence that soft gamma repeaters (SGRs) are neutron stars which 
experience frequent starquakes, possibly driven by an evolving, ultra-strong 
magnetic field. The empirical power-law distribution of SGR burst energies, 
analogous to the Gutenberg-Richter law for earthquakes, exhibits a turn-over 
at high energies consistent with a global limit on the crust fracture size. With 
such large starquakes occurring, the significant excitation of global seismic 
oscillations (GSOs) seems likely. Moreover, GSOs may be self-exciting in a 
stellar crust that is strained by many, randomly-oriented stresses. We explain 
why low-order toroidal modes, which preserve the shape of the star and have 
observable frequencies as low as ~ 30 Hz, may be especially susceptible to 
excitation. We estimate the eigenfrequencies as a function of stellar mass and 
radius, and their magnetic and rotational shiftings/splittings. We also describes 
ways in which these modes might be detected and damped. There is marginal 
evidence for 23 ms oscillations in the hard initial pulse of the 1979 March 5th 
event. This could be due to the 3^0 mode in a neutron star with B ~ 10 14 G or 
less; or it could be the fundamental toroidal mode if the field in the deep crust 
of SGR 0526-66 is ~ 4 x 10 15 G, in agreement with other evidence. If confirmed, 
GSOs would give corroborating evidence for crust-fracturing magnetic fields in 
SGRs: B £ 10 14 G. 

Subject headings: stars: magnetic — stars: neutron — stars: oscillation- 
X-rays: stars 
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1. Introduction 

Global shear modes in the Earth were first detected following the devastating Chilean 
earthquake of May 1960 (Benioff, Press & Smith 1961). The fundamental toroidal mode, 
which for m = involves a simple, alternating, periodic twisting of whole hemispheres 
in opposite directions, has a period of 43 minutes in the Earth. In general, the toroidal 
modes are pure shear deformations, divergence- free, with no radial components. They 
are referred to by the notation ct n , where the overtone number n is the number of radial 
nodes in the eigenf unctions (e.g. McDermott, Van Horn & Hansen 1988; hereafter MVHH). 
The fundamental mode is 2^0, because £ = 1 oscillations would violate angular momentum 
conservation. Toroidal modes up to £ = 60 have been detected in the Earth, with 
P(mto) = 2.3 minutes (Lay & Wallace 1995). 

Any condensed object with a solid component can support shear modes. For example, 
old white dwarfs with crystalline cores have -P^o) ~ 12 s (Hansen & Van Horn 1979), but 
in this case, the solid core lies within an opaque, liquid mantle, so detection is problematic. 
Neutron stars, which have fractionally-thin solid crusts floating on mantles of superfluid, 
offer better prospects for supporting observable shear modes. Because the crusts are thin 
and subject to little friction from below, large-amplitude shear oscillations could be excited 
with much less energy than would be required in true solid bodies. Moreover, substantial 
modal deformation occurs at, or very near, the observable surface. 

The toroidal (or torsion) modes, which preserve the star's shape, may be especially 
easy to excite via starquakes, because the restoring force for these modes is entirely due 
to the relatively weak Coulomb forces of the crustal ions as given by the shear modulus: 
li ~ (Ze) 2 /a 4 , where a is the ion spacing. Other internal modes (e.g., spheroidal modes in 
the crust) involve bulk compression and vertical motion, which have to do work against 
the much stronger degeneracy pressure of electrons in the outer crust and free neutrons in 
deeper layers, as well as gravity. This implies: (1) toroidal shear deformations require much 
less energy than do radial or compressional deformations of comparable amplitude; and (2) 
low-order t-modes also have unusually low frequencies. Since the damping rate (e.g. via 
coupling to Alfven modes in the magnetosphere or other mechanisms; §3 below) generally 
increases steeply with oscillation frequency, both these circumstances favor the significant 
excitation of low-order t-modes. 

For many — or all — of these reasons, shortly after the discovery of pulsars, Ruderman 
(1968) proposed that the radio pulsations found by Hewish et al. (1968) were due to the 
fundamental toroidal mode, 2 to, m a neutron star. We now know that spindown-driven 
seismic activity is possible in ordinary pulsars (e.g., Ruderman 1991), but it has not proven 
sufficient to excite detectable crust vibrations. However, there is another class of stars 
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in which Ruderman's prescient proposal might more readily apply: soft gamma repeaters 
(SGRs). 

A variety of evidence seems to indicate that SGRs are neutron stars which experience 
strong and frequent starquakes (e.g., Thompson & Duncan 1995, 1996; hereafter TD95 
and TD96). Comparative analyses of SGR and earthquake records (Cheng et al. 1996) 
indicate that the relative distribution of SGR burst energies and waiting times, and their 
correlations, are consistent with crust fracture/slippage events, but not with any known 
accretion- or nuclear-powered phenomenon. This accords with the hypothesis that SGRs 
are neutron stars with ultra-strong magnetic fields, B 3> m^c 3 /eh = 4.4 x 10 13 G. Under 
this hypothesis, SGR bursts are the observable signature of starquakes driven by the 
stresses of an evolving magnetic field; and the exceptionally energetic 1979 March 5 event is 
identified as a magnetic flare (Duncan & Thompson 1992, hereafter DT92; Paczyhski 1992; 
TD95). Such magnetically-active neutron stars, or magnetars, could form as a result of a-Q 
dynamos in hot, newborn neutron stars undergoing convective mixing (e.g., Keil, Janka & 
Miiller 1996) if they are born rotating more quickly than the relevant "dynamo number" 
threshold (DT92; TD93; DT96). It is even possible that some magnetars form without 
strong convection (e.g., following rotation-supported core bounce) because differential 
rotation, which is generic for high-angular momentum proto-neutron stars, could alone 
generate fields ~ 10 17 (P rot /l ms) _1 G, where P rot is the rotation period (DT92; TD93; 
DT96). Seven different estimates of the field in SGR 0526-66 indicate B £ 10 14 G (§1 in 
TD95). 

We now consider the physics of global seismicQ oscillations (GSOs) in SGRs. We will 
focus on t-modes, although our results on excitation and damping could apply to other 
crustal modes (e.g., s-modes) as well. 

2. Eigenfrequencies and Strain Amplitudes 

The fundamental toroidal mode's period is roughly given by the time for a shear wave 
to travel around the star: P ~ R/V^, where = (/i/p) 1 / 2 is the shear wave velocity. 
This crude estimate ignores vertical structure in the crust; but this is a surprisingly good 
approximation because varies only weakly with depth (Ruderman 1968; Strohmayer 
et al. 1991). Accurate mode calculations have been done by Hansen & Cioffi (1980) 



^^Here, seismic refers to waves excited in the crust by starquakes, in analogy to the common geophysical 
usage. This should not be confused with seismological, which could refer to any internal mode probing the 
interior. 
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and MVHH. Separation of variables for the shear wave equation, subject to the t-mode 
conditions u r = and V ■ u = 0, gives solutions of the form 

W{r) dY em dY em 

where {u r ,ue,u^) are displacement amplitudes in spherical coordinates, and Y£ m (6,<j)) are 
spherical harmonics. For n = modes, the radial eigenf unctions W(r) vary little with 
depth in the crust (MVHH). 

Results of MVHH's model calculations are given in the first and third columns of Table 
1. Note that the neutron star radii listed here are smaller than generally thought realistic 
for their respective masses, because MVHH used an interior (sub-crust) equation of state 
that is probably unrealistically soft. We now seek to estimate of the observable eigenperiod 
for a canonical M = 1.4 M , R — 10 km neutron star (center column of Table 1), as well 
as its dependence on M and R. To do this, note that the structure of the crust depends 
upon M and R almost entirely through the action of the surface acceleration of gravity; 
hence so will the effective (depth-averaged) V^. The quantity P/R oc V~ x (fourth row of 
Table 1) should therefore be a function of g — GM/R 2 (fifth row, tabulated in units of 10 14 
cm s -1 ); so we interpolate over g to we estimate P/R and P (center column). Two further 
corrections are needed. Improved calculations of // (Strohmayer et al. 1991), including 
directional averaging (Ogata & Ichimaru 1990), imply that P oc /i -1 / 2 is actually larger 
by 1.24 than found by MVHH, as quoted in Table 1. Gravitational time dilation further 
increases the observable period. We find: 

x r. 0.87 + 0.13M 14J Rr 2 

PUt ) = 33.6 R w - 1° ms, 2 

v ; (1.71 -0.71ML4.Rro) 1/2 

where R w = (R/10 km) and M 1A = (M/1.4M ). 

For higher-order modes, the period goes as P(eto) oc [£(£+ l)] -1 ^ 2 (e.g., Hansen & Cioffi 
1980), so the predicted mode spectrum, in milliseconds, is roughly P pa 34, 24, 18, 15 . . . 
Each of these modes is (2£ + l)-fold degenerate in m. (In §4 we consider magnetic and 
rotational shifting and splitting of the eigenfrequencies.) Overtone modes (n > 0) have 
much shorter periods, of order a few milliseconds or less (e.g., MVHH). 

The mode energy scales roughly as E m oc R 4 /M since the elastic deformation energy 
per unit area of crust goes as ~ g^ 1 in hydrostatic equilibrium [eq. (6) below]. Thus the 
quantity (EM/R 4 ) is nearly constant in MVHH's models (bottom row of Table 1), where 
E = E m [W(R)/R}~ 2 (sixth row; tabulated in units of 10 48 ergs). Interpolating over g 
again, and including the correction factor to /i, we find 



E m = 5 x 10 



47 



W(R) 
R 



2 



i?io M^ 1 (0.77 + 0.23M L4 i?r 2 ) erg (3) 
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for the 2 to mode energy. The surface displacement amplitude for m = is then given by 
eq. (1): 

U<t> = 440 E\{ 2 i?" 1 M\[l sin 29 cm, (4) 

where Ku = (i? m /10 41 erg). [For simplicity, we omit the fine adjustment factor 
(0.77 + 0.23Mi. 4 i?7 2 ) -1 / 2 in this and the following equation.] The amplitude of the modal 
strain is 

i) m = 0.9 x 10~ 3 EH 2 Rio M\[l sin 2 6. (5) 
For m/0, the angular distribution of strain is different but the amplitude is comparable. 

Equation (5) is a remarkable result, because: (1) the observed, prompt X-ray energy 
a single bright SGR burst is ~ 10 41 ergs (e.g., Norris et al. 1991) so the energy of a major 
starquake must be comparable or larger than this; and (2) the static critical (yield) strain 
of a neutron star crust is probably in the range ip cr ~ 10~ 3 , and perhaps less (Smoluchowski 
& Welch 1970; Ruderman 1991). [Although if) cr ~ 10~ 2 in a pure bcc Coulomb crystal, 
realistic crustal material is almost certainly weakened by lattice dislocations.] Thus if a 
seismic episode, involving multiple starquakes, imparts to the mode even a minor fraction 
of the X-ray energy of a single bright burst, then the modal strains themselves will 
be large enough to trigger starquakes at other sites in the crust that have pre-existing 
(magnetically-induced) strains significantly less than ip cr . These triggered events can give 
impetus to the mode, as we now explain. 



3. Mode Excitation, Damping and Detectability 

The energy of an event in which critical strains are relieved throughout an area A of 
the crust is £ = (1/2) 4> 2 r At; J fidz, where the integral runs over the depth of the fracture, 
and £ > 1 accounts for released magnetic (as opposed to elastic) energy. (This satisfies 
£ < (Airfx/B 2 ) ~ 10-Bf 5 2 except when there is significant reconnection; §2.2 in TD95.) Using 
the equation of state of Negle & Vautherin (1973) to evaluate the integral down to the base 
of the crust, we find 

£ = 4.2 x 10 38 l\ tf r3 i Rj Mr} erg, (6) 

where ^3 = (VW 10~ 3 ) and we take A = lj for a fracture of length If = 1 Zi km. This 
result scales as g^ 1 oc i? 2 M _1 because dz = (pg)~ 1 dP in hydrostatic equilibrium. 

The empirical burst energy distribution for SGR 1806-20 is a power law with index 
7 = 1.66 (Cheng et al. 1995), analogous to the Gutenberg- Richter law for earthquakes 
(Gutenberg & Richter 1956). This is a natural consequence of self-organized criticality (e.g., 
Bak 1996). Unlike the case of earthquakes, the SGR burst distribution turns over at high 
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S, with an apparent upper limit 8 max = 3 x 10 41 Df 4 ergs, for SGR 1806-20 at a distance 
D = 14D 14 kpc from Earth (Corbel et al. 1997). If this is due to a maximum fracture scale, 
lf,maxi then 

(l f>max /2nR) = 0.4 r 1/2 r 1 ' 2 Rio Ml 1 ] Du , (7) 

where / <^ 1 is the fraction of the released energy that powers burst emissions. This 
suggests that the observed turn-over in the SGR burst energy distribution is due to global 
limits on the fracture size. For comparison, the largest earthquakes yet measured by modern 
seismographs have dimensions (lf/2irR <s ) < 10~ 2 (Scholz 1990). 

Global fractures propagate across the crust on time scales ~ R/V^, comparable to 
the period of low-order seismic modes, which are thus readily excited. Once GSOs begin, 
new fractures can be triggered, but always as the modal strain is locally increasing with a 
component in the direction of the pre-existing (magnetic) strain. Each of these triggered 
fractures can give an energy boost to the mode, for reasons that may be understood via 
a simple 1-D analogy. Imagine that a rubber band is stretched to near its breaking point 
and held. This requires energy, say from person A who stretches the rubber band; thus 
A plays a role analogous to the magnetic field. A second person, B, grasps the band just 
within A's fingers and imparts a further, small-amplitude periodic stretching, as done by 
the GSO. When the band breaks, it snaps against B's fingers, giving an impulse in the 
direction of B's motion (a mode-boosting impulse). The energy released is much greater 
than the energy of stretching by B (the mode); it ultimately comes from A (the magnetic 
field) as stored in the elastic band (the crustal strain). After the band breaks, person A 
might also do additional work in dragging B's fingers apart. This is the effect quantified by 
£ above. It depends upon the fault slippage, which is difficult to estimate (§2.2 in TD95). 

This analogy suggests that GSOs are potentially self-exciting in a stellar crust with 
many randomly- oriented strains. Indeed, for a power-law distribution of pre-existing strains, 
[d 2 N jdl dip] oc l~ a if}®, and for an excitation efficiency e = e Q l\ , self-sustaining mode growth 
occurs to modal strain amplitudes ipm ^ (<£m • 4>cr) if e G > XV 9 ™ (1 + y?m) _1 (1 — <£>m) _ ^ Q+1 \ 
where 

X = 0.47 ■ (27)"' ( { J^ { ~° a) ) ^3 A" 1 < +2 M^ /2 , (8) 

and where A is the fraction of the star's surface area that is strained. The SGR 
"Gutenberg-Richter Law" (Cheng et al. 1995) implies a = 2.32 if most fractures occur 
when ip m = ip cr — ip <C ip cr . As long as (5 > a — 3 ~ —0.68, large fractures dominate 
the excitation. This is probably true since large fractures tend to be more efficient at 
excitation than small, i.e., (5 > 0. For (3 = 1 = a and ip m <^ 1, the efficiency criterion is 
e D ^ 0.002 (</? m /0.1) ip cr3 A -1 £ -1 ^ 2 Rf M^ 2 . This suggests that terminal mode amplitudes 
~ ip cr could be attained by self-sustained mode growth. 
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When excitation is dominated by large fractures, it is possible that only a few low-order 
modes are substantially excited. For high-(£, m) modes, a randomly-oriented fault line can 
traverse modal zones with changing senses of strain, which tends to hamper excitation. 
As starquakes are triggered, the mode axis will wander (since excitation amplitudes add 
in phase like axial vectors) allowing strains to be released all over the crust. Thus modes 
which happen to be excited by a spontaneous starquake at the start of a seismic episode 
could grow rapidly, leaving little elastic/magnetic energy to be tapped by other modes. 

How might GSOs be detected? SGR burst emissions could be significantly modulated 
at the mode period. For example, a large crust fracture, once began, could extend in 
length each time the mode stresses it the right way, producing new pulses of Alfven 
waves, pair creation, and observable X-rays during each during each GSO cycle (cf. §2 
in TD95). Alternatively, there may be detectable discreteness in the time-spacings of 
the aftershock-like events following a bright burst (Kouveliotou et al. 1996), if these are 
triggered by the mode rather than by the general redistribution of stress. Since triggering 
usually occurs near but just before the two instants of maximum deformation in the mode 
cycle, the expected time intervals between mode-triggered bursts are At ~ i](P/2), where rj 
is an integer. This is especially true for At ^ P. When At ^> P, especially over intervals 
containing bursts, modal phase drift would mask the effect. 

How do GSOs damp? There are several possible mechanisms; here we will concentrate 
on fault line reslippage and refracturing. Like the Earth, a neutron star is probably laced 
with fault lines, along which the effective critical strain is smaller than usual: ipcrj < ipcr- 
These may be sites of recent fractures that have not fully healed (in which case ip cr j is 
actually the threshold for slipping against friction rather than a true yield strain; Sholtz 
1990), or areas where magnetic stresses have a history of accumulating and weakening 
the crust. Whenever the modal strain amplitude along the fault line, ip m j, exceeds ip C rj, 
then slippage occurs, as often as ~ 2ip m j/tp cr j times per cycle when this number exceeds 
two. This "sticky" motion reprocesses mode energy, mostly into seismic waves of peak 
frequency ~ V^/Z s , which is in the ~ kilohertz range for l s ~ 1 km, comparable to the 
depth of the crust. Such high-frequency waves couple strongly to Alfven modes in the lower 
magnetosphere (TD95; Thompson & Blaes 1997), while low-frequency ones are internally 
reflected (Blaes et al. 1989; paper in preparation). Each reslippage/refracture event thus 
injects energy into the magnetosphere, maintaining an X-ray emitting plasma there. The 
resultant mode damping time t s , is difficult to estimate, but if t s < 130 F 7 -1 (_D/8kpc)~ 2 
s, where F th = 10~ 7 F 7 erg cm" 2 s' 1 is the X-ray detector threshold, then the energy of 
GSO damping would be detectable as quasi-steady X-rays during active episodes, fading 
after bursts cease. Such emissions could be detectable even if many different seismic modes 
are excited. 
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4. Magnetic Frequency Shifts and the March 5th Event 

If a tangled magnetic field of r.m.s. strength B is embedded throughout in the crust, 
then magnetic tension augments p, and 

P^P [l + (B/B,f]- l /\ (9) 

where P Q is the non-magnetic period (eq. [2], in the case of 2^0), and 

B^ = (471-yu) 1 / 2 a: 4 x 10 15 p54 4 G, based on a power-law fit to the deep-crust 

(p > Pdrip = 4.6 x 10 11 gm cm~ 3 ) equilibrium composition found by Negle & Vautherin 

(1973). Note that the mode energy (eq. 3) also shifts, upward by a factor [1 + (B/B^) 2 ]. 

Even if B <sC B^ in the deep crust, near the star's surface, at densities p < 4 x 10 10 B\f 4 
gm cm -3 for B u < 4, magnetic restoring forces exceed elastic ones. This perturbs the 
eigenfunctions in the outer layers of the crust, but it has little effect on the observable 
GSO period since the great preponderance of mode deformation energy (or kinetic energy) 
resides below this (Carroll et al. 1986). Thus it is appropriate to take pu ~ 1 in the formula 
for B^ when evaluating eq. (9) [except in the Alfven limit, B ^> B^, for which eq. (9) no 
longer applies]. 

Magnetic frequency shifts also can split the eigenfrequencies, resolving the (2£+ l)-fold 
degeneracy in m. Equation (9) is precisely accurate only in the case of a disordered 
field uniformly penetrating the star's crust. A realistic field with global structure will 
produce m-dependent period shifts of order (AP/P) ~ [B/B^) 2 for B <C B^, where the 
precise splittings depend upon field geometry and mode orientation. Rotational splittings 
AP/P = [m/£(£+ 1)] (P/Prot) also occur (Lapwood & Usami 1981; Strohmayer 1991). 
However, for 2 t , this is only AP/P = 7 x 1CT 4 m (P rot /8.0 s)' 1 , where we scale to the 
rotation period of SGR 0526-66 as inferred from the 1979 March 5th event. For comparison, 
rotational splittings of t-modes are ~ 10 times larger in the Earth, yet they have not been 
detected (Lay & Wallace 1995). 

The exceptionally energetic 1979 March 5th event was probably due to a large-scale 
instability in a magnetar (§2.1 of TD95). A likely candidate is the Flowers- Ruder man 
instability, which involves a global crust fracture and a (small-angle) relative twist of two 
sections of crust, diminishing the exterior dipole moment as reconnection powers the flare 
(Flowers & Ruderman 1977; §14.2 and §15.2 in TD93). This would excite very strong GSOs. 
Indeed, there is evidence for a 23-ms periodicity in the hard initial pulse of the event (Barat 
et al. 1983). This could be due to 3 t oscillations in a neutron star with B ^ 10 14 G; or it 
could be due to the fundamental, 2 to, if the field in the deep crust is B = 1.06 B^ ps 4 x 10 15 
G, for a 1.4 M , 10 km neutron star (eqs. [1] and [9]). Such strong fields in the deep crust 
of SGR 0526-66 have previously been inferred for other reasons (TD95; TD96). Note that 
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several modes would probably be excited by such a catastrophe, and the ongoing release of 
energy would also affect crustal motions, so a complex (non-sinusoidal) light curve is to be 
expected. 

In conclusion, if SGR bursts are due to starquakes, then strong seismic waves almost 
certainly vibrate SGR crusts during active episodes. Whether a few low-order GSOs can 
attain sufficient amplitudes to have detectable consequences is less certain. We have given 
three reasons why this is at least possible: (I) global-sized fractures (eq. [7]) naturally excite 
low-order modes; (2) mode damping, via Alfven wave loss and other mechanisms, declines 
steeply with frequency; and (3) the "self-exciting" effect could drive runaway growth of 
single modes. 

If the presense of GSOs in SGRs, as suggested by the 1979 March 5th data, is verified 
by future observations, this will add to the (already considerable) evidence for seismic 
activity in these stars. Magnetism is a natural candidate for driving this activity, especially 
in light of the other evidence for very strong fields in SGRs (as summarized in §1 of TD95). 
Rotational energy may trigger some starquakes and glitches in ordinary pulsars (e.g., 
Ruderman 1991), but this could not power SGR activity in a neutron star with an 8.0 s 
period like SGR 0526-66. Magnetars seem necessary, because typical pulsar-strength fields 
could not drive strains in excess of tp cr and fracture the crust. Hooke's Law requires 



where p i4 ~ 1 because only deep fractures have sufficient energy to power bright SGR 
bursts (eqs. [6]-[7]). 
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Table 1. Eigenperiod and Energy of the Fundamental 
Toroidal Mode 



M/M Q 0.503 1.4 1.326 

R (km) 9.785 10 7.853 

P (ms) 18.54 20.7 17.32 

P/R (ms/km) 1.89 2.07 2.21 

gu 0.701 1.87 2.87 

E i8 1.60 0.73 0.331 

E 48 M 1A /Rf 0.63 0.73 0.82 



Note. - Columns 1 and 3 show results for models 
NS05T8 and NS13T8 of McDermott, Van Horn & Hansen 
(1988), which employ a soft interior equation of state. 
Column 2 gives extrapolated values for a canonical model 
star. Further corrections are needed before comparing with 
observations, as explained in the text. 



